Specular Andreev Reflection in Inversion-symmetric Weyl-semimetals 
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The electron-hole conversion at the normal-metal superconductor interface in inversion-symmetric 
Weyl semimetals is investigated with an effective two-band model. We find that the specular Andreev 
reflection of Weyl fermions has two unusual features. The Andreev conductance for s-wave BCS 
pairing states is anisotropic, depending on the angle between the line connecting a pair of Weyl 
points and the normal of the junction, due to opposite chirality carried by the paired electrons. For 
the Fulde-Ferrell-Larkin-Ovchinnikov pairing states, the Andreev reflection spectrum is isotropic 
and is independent of the finite momentum of the Cooper pairs. 

PACS numbers: 74.45. +c, 73.23.-b 



In recent years, a great progress in condensed matter 
physics has been made by the discovery of topological 
insulators [IJ, which have a bulk energy gap and gap- 
less surface states protected by the time reversal sym- 
metry. Such a phase is characterized by the topological 
invariants rather than any order parameters. Currently, 
the topological matter is further extended to the Weyl 
semimetals (WSMs) [2], where the bulk energy band is 
gapless and the Weyl points (WPs) emerge separately in 
the momentum space in pairs carrying opposite chiral- 
ity. Electrons around the WPs can be well described by 
the two-component Weyl equations. In order to realize 
this fascinating phase, either the time-reversal symmetry 
[2H5] or the inversion symmetry [6] needs to be broken. 
The WSM phase is protected by the band topology and 
therefore is robust. The Dirac type dispersion cannot be 
gapped unless there exists coupling between the pair of 
WPs of opposite chirality [2j [3] . 

When the superconducting pair potential is introduced 
into those topological matters through proximity effects 
[7] or intrinsic phonon-mediated attractive interaction 
[8], a superconducting phase with non-trivial topology 
can emerge. Taking the topological insulator as an ex- 
ample, the s-wave pair potential can result in a topologi- 
cal superconductor with Majorana surface states [HO [5]. 
Recently, such a route has been adopted by Meng et al. 
[§] and Cho et al. [TU] in the study of superconductivity 
in inversion-symmetric WSMs. They found a novel BCS 
pairing state which is an electronic analogue of 3 He-A 
phase and the Majorana modes can emerge under certain 
conditions [HI [TU]. Apart from the BCS pairing states, 
the separation of WPs in momentum space also provides 
an opportunity for the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) pairing states [H][T2], where the Cooper pairs 
carry a finite momentum. In contrast to the conventional 
case, the FFLO state in WSMs is robust against the weak 
disorder [TU] . 

The study on the superconducting WSMs is still in the 
preliminary stage and one central question to be resolved 
is how to detect them. One of the effective methods 
for detecting pairing potentials is through the Andreev 
reflection (AR) [13j at the normal-metal superconductor 



FIG. 1. (Color) Illustration of the NS junction composed 
of WSMs. The black and white balls represent the elec- 
trons and holes, respectively. An electron of incident an- 
gle (9, <p) from the normal-metal side (the left) is scattered 
into four branches: the normal reflection, the specular AR, 
the electron-like quasiparticle transmission, and the hole-like 
quasiparticle transmission, with amplitudes r e , r^, t e , and th, 
respectively. Line P-P+ makes an angle a with the normal 
of the interface. 



(NS) interface. The AR spectrum has been successfully 
used to probe both the conventional s-wave [J_D and the 
unconventional anisotropic [T5HT7] superconductors. The 
recent progress of the studies on the AR also reveals some 
unusual features, such as the specular AR in the graphene 
[T5] and the resonant AR induced by the Majorana modes 

ping. 

In this paper, we employ an effective two-band model 
to investigate the specular AR at the NS interface in 
inversion-symmetric WSMs with either the BCS or the 
FFLO pairing states. It is found that the normalized 
subgap conductance for the s-wave BCS pairing states is 
a constant of 1.726, and the effective gap is anisotropic, 
which is a sinusoidal function of the angle between the 
axis connecting two WPs and the normal of the junction. 
For the FFLO pairing states, the AR spectrum exhibits 
an unexpected isotropic behavior, being nothing to do 
with the finite momentum of the Cooper pairs. 

The NS junction under consideration is sketched in Fig. 
[TJ The normal of the junction is along the z-axis and the 
interface is located at z = 0. The superconducting WSM 



can be prepared by doping [TO], which drops the WPs 
well below the Fermi level. The doping level should be 
chosen properly, so that the electrons at the Fermi level 
have a large density of states for the phonon-mediatcd 
condensation but they can still be described by the Weyl 
equation. 

The minimal model for inversion-symmetric WSMs is 
described by an effective two-band Hamiltonian (in units 

of n = l) E3 

a=± (1) 

h±(q,<r) = v (q x a x + q y a y T^z^z) - U<d{z), 

where c a — {c a ^,c a i) T is the electron operator with 
a = ± denoting two Weyl nodes, respectively, vq is the 
Fermi velocity, q is the momentum measured from the 
WPs, a is the Pauli matrix of spins, and 6(2) is the step 
function. Each Weyl node contains two cylindrical three- 
dimensional cones that touch at the WP. In this model, 
a pair of WPs are located at two inversion-symmetric 
points, .P+ = P and P = P. In general, the axis 
connecting two WPs, line P-P+, can make an angle of 
a with the normal of the junction, as shown in Fig. [TJ 

cos a — sin a\ 



Therefore, the rotation operator R x = 

\sma cos a 

is incorporated in Eq. (nj), with which a vector is ro- 
tated around the y-axis by an angle a. When a = 0, 
line P-P+ is along the z-axis and the Hamiltonian in 
Ref. [TU] is recovered. The energy in Eq. ([!]) is measured 
from the Fermi level. In the normal-metal side (2 < 0), 
there is no band shift and the WPs are at the Fermi level. 
In the superconductor side (z > 0), there is a large band 
shift U caused by doping, so that the WPs are well below 
the Fermi level. Although the effective Hamiltonian Eq. 
(fTl) contains only one pair of WPs, it can well describe 
the low energy excitations of two proposed candidates for 
WSMs based on topological insulators [H [5] . 

The electrons in WSMs bear the spin-momentum lock- 
ing features and nodes P± carry opposite chirality. As 
a result, there are two types of pairing mechanisms [TU] : 
one is the inter-node BCS pairing for which the paired 
electrons are from two Weyl nodes of opposite chirality 
and the momentum of the Cooper pair is zero, the other 
is the intra-node FFLO pairing for which the paired elec- 
trons are from only one Weyl node and the momentum 
of the Cooper pair is finite (2P± ) . The Hamiltonian for 
both pairing states can be expressed as 



i? p B air = J2 A ( z ) c Uv) c -ai(-q) + h.c. 



q , a 



^ = J2 A ( z ) c U<i) c ii(-i) + h - c -> 



(2) 



where superscripts B and F correspond to the BCS and 
FFLO pairing states, respectively, and A(z) = A6(z) is 
the s-wave pairing potential, which is chosen to be real. 



After the diagonalization of the Hamiltonian H = Hq- 
•Hpajjf j the energy excitations are given by 



eB = y{vo\q\ ± VU 2 + A 2 cos 2 0) + A 2 sin 2 /?, 

E F = ^(v \q\±U) 2 + A 2 , 



(3) 



where /3 is the azimuthal angle of the momentum relative 
to line P-P + . Due to the different manners of spin- 
momentum locking, the electrons from two Weyl nodes 
cannot fully paired [10] , The excitation spectrum E B 
is consequently characterized by an anisotropic effective 
gap I A sin j3 | . The nodes of E B exactly lie in line P-P+ 
(/3 = 0,7r). On the contrary, the Cooper pairs in FFLO 
superconductors are formed within the same Weyl node 
and the excitation spectrum E F is fully gapped. 

The Bogoliubov-de Gennes (BdG) Hamiltonian in the 
real space can be obtained by performing a Fourier trans- 
formation on Eqs. pj and (pi). In the Nambu representa- 



tion where tp 



bfo(r),ipi(r),il>l(r), 



V>j(r)] T , the BdG 



Hamiltonians for both pairing cases are given by 



-"BdG 



U F — 
-"BdG — 



h' ± (-iV-P ± ,cr) 

AO(z) 

fc' ± (-iV-P ± ,«7) 



AG(z) 
-h'(-iV-P±,a) 



A, 



,2iP+ 



-9(z) 



Ae- 2iP± - r Q(z) -ti ± (-iV - P T ,a) 



(4) 



where h'±(V,cr) = /i±(i? _1 V, Br l cr). The pairing po- 
tential exhibits a spatial modulation for the FFLO su- 
perconductors and is a constant for the BCS ones. 

Considering an electron of energy E and of incident 
angle {9, tp) coming from the normal-metal side, it can 
be normally reflected as an electron or Andreev reflected 
as a hole, or transmitted into the superconductor side 
as the electron-like or hole-like quasiparticles, as shown 
in Fig. [T] As the Fermi level is located exactly at the 
WPs in the normal-metal side, the incident electron and 
the reflected hole must belong to the conduction and the 
valence band, respectively, resulting in the specular AR, 
which is similar to what occurs in the graphene [T5] . 

For the BCS pairing WSM, an incident electron in node 
P + can be converted into a hole in node P_ via the AR 
process. The momenta of the electron and hole are both 
around P + , so that the BdG equation can be solved with 
the ansatz tp = (fe lP+ ' r ,ge lP +' r ) T , where the / and g 
are the electron and hole components of the wave func- 
tion. For the FFLO pairing WSM, an incident electron 
in node P_ is converted into a hole still in node P_ . The 



momenta of the electron and hole are around P_ and P 



+ • 



respectively, and the BdG equation can be solved with 
the ansatz V = {fe lP - r ,ge lP + r ) T . 

For the BCS pairing case, the wave function in the 
normal-metal side ^^ and that in the superconductor 



side Wg are given by 






b ie lqcosez + rfb 2 e- iqcosez + r^b 3 e~ iq 



cos 6z 



*|=tf6 4 e^+ z +^& 5 e I 



(5) 



where h = [s(a, 9), s(a - tt/2, (9), 0,0] T , 6 2 = &i(tt - 



[0,0,e-^/ 2 cos(9/2, 



Jip/2 



sin(9/2] T , 6 4 



[sin2a, -2cos 2 a, Asin2a/(£' + E~/ B ) 7 0] T , and b 5 = 
[sin 2a, 2 sin a, 0,2Asin 2 a/(E — _E7s)] T are the basis 
functions for the incident electron, the normal reflected 
electron, the Andreev reflected hole, and the electron- 
like and hole-like transmitted quasiparticles, respectively, 
with s(a,9) = cosasin(6»/2)e" lv / 2 - sinacos(6»/2)e 1 ^ 2 , 
and 7s = \J E 2 — A 2 sin 2 a/E. The scattering coeffi- 
cients for each process are denoted by rf h and tf h , re- 
spectively. The wave vectors are obtained as q — E/vq 
in the normal-metal side, and qf ± = (±C/ + E^b)/vq for 
the electron-like and hole- like quasiparticles in the super- 
conductor side. The rapid oscillation factor e lP +' r in the 
ansatz and the transverse plane wave e I,? ii' r are omitted 
in Eq. (|5|) for simplicity, as the transverse momentum q\\ 
is conserved during the scattering. In Eq. ([5]) , the heavy 
doping condition U ^ E, A is also utilized. As a result, 
there is a large mismatch of wave vectors between two 
sides, and the quasiparticles in the superconductor side 
are nearly perpendicularly transmitted. 

For the FFLO pairing case, the wave function can be 
obtained in a parallel manner, and reads 



*£ = /ie l9COs92 +rf/ 2e -^ c 



F i 

r h h& 



-iq cos Oz 



(6) 



where the basis functions are written as /1 = 
[( e -^/ 2 cos6»/2,e^/ 2 sin6l/2)e lP -- T ',0,0] T , / 2 = /i(tt - 
9), h = [0,0, (e-^/ 2 cos6»/2,-e^/ 2 sin6»/2)e 4P +- r ] T , 
h = [e iP - r ,0,e ip +- r E(l - 7 f)/A,0] t , and f 5 = 
[0,e lP - r ,0,e lP + r E{l + 7f)/A] t , respectively, with 
jp = y/E 2 — A 2 /E and the wave vectors qf ± = (±U + 
Ejf)/vq. One finds that the rapid oscillation factors 
e iP±-r j? Qr ^ e e i ec t ron anc | hole components of the wave 
functions are different for the FFLO superconductor. 



By matching the boundary condition ^ 



B(F) 

N 



= <J 



B(F) 



at the interface z ~ 0, one obtains the scattering ampli- 
tudes as 

(Asina/.E)cos# 



■f = : 

cos </?(7_b + cos ff) — i sin ip(l + 73 cos 9) ' 

sin 6(— jb cos ip + i sin tp) 



„B 



B 



COS ip(j 

(A/E) cos 9 
1 + 7F cos 9 ' 
sin# 



cos( 



i sin ^(1 + jb cos( 



(7) 



e 1 + 7f cos 9 

With the help of the Blonder-Tinkham-Klapwijk 
(BTK) formula [14] . for a specific incident angle, the 
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FIG. 2. (Color) Contour plot of the normalized differential 
conductance of the NS junction with the BCS pairing states, 
as a function of eV and a. 



differential conductance is determined by the BTK co- 
efficient 



S*(E, 9,<p) = l+ K(E, 9, p)\ 2 - \rl(E, 9, <p)f 



(8) 



with i = (B, F). In units of the conductance of the junc- 
tion in the normal state (A = 0), the normalized conduc- 
tance is given by 



r-TT/2 



G l (eV) 



ft'* d9 sin 29 J dtp&ieV, 9, tp) 



r ir/2 



,-2jr 



/*'' d9 sin 29 /„ dp5*(eV, 9, p)\ A=0 



(9) 



where eV is the bias voltage. 

Firstly, we consider the AR spectrum for the BCS pair- 
ing states. Although the pair potential of an s-wave BCS 
pairing WSM is a constant, the electrons with different 
spin-momentum locking features can only partly paired 
according to their momentum direction, leading to an 
anisotropic effective gap |Asin/3| in energy excitation 
E B . Such an anisotropic gap is also exhibited in the 
AR spectrum. Under the heavy doping condition, the 
quasiparticles are nearly perpendicularly propagated in 
the superconducting WSM. The azimuthal angle of the 
momentum relative to line P-P+ can be approximated 
as a. Therefore, one finds that the AR coefficient rf L is 
proportional to A sin a and the jb is also determined by 
A sin a. When a = or 7r, the effective gap is zero and 
the AR is totally suppressed. 

The dependence of the normalized conductance on 
both the bias voltage and angle a is plotted in Fig. [2] A 
notable characteristic in Fig. [2] is the yellow region of a 
contour of |A sina|, in which the normalized conductance 
is a constant of 1.726. Such a constant behavior can be 
understood as follows. When the incident electron is of 
an energy below the effective gap, the transmitted wave is 
evanescent and jb is pure imaginary, so that the subgap 
conductance is determined by the BTK coefficient 



S° = 2|r 



B|2 

h I 



2 cos 2 



is 2 9 cos 2 (ip — if) 4- sin (tp — r/) ' 



(10) 
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FIG. 3. (Color) Plots of the normalized differential conduc- 
tance as a function of eV. The solid line is corresponding to 
the FFLO pairing states. The dashed, dotted and dot-dashed 
lines are corresponding to the BCS pairing states with differ- 
ent a, respectively. 



with r\ = arctan(v A 2 sin 2 a — E 2 /E). For the per- 
pendicular incidence (9 — 0), the normal reflection of 
Weyl fermions is completely suppressed due to the spin- 
momentum locking, and BTK coefficient S< is a constant 
of 2, as in the conventional NS junction in the transpar- 
ent limit [14] . For the oblique incidence, there is a finite 
normal reflection caused by the large momentum mis- 



match at the interface, and 5v in Eq. (jlOh is less than 2 



and depends on the energy in general. Interestingly, the 
energy factor in S< is always bound with the incident 
angle (p. Therefore, after an integration over incident an- 
gles and normalization, the subgap conductance becomes 
a constant of 1.726, independent of energy. 

Secondly, we consider the AR spectrum for the FFLO 
pairing states, which is plotted in Fig. [3] as the solid 
line. Although the Cooper pair in an FFLO pairing WSM 
bears a finite momentum of 2P± and the pair potential is 
modulated in the spatial space, the AR spectrum is still 
isotropic and has nothing to do with P± . This is due to 
three reasons: a) In contrast to the conventional FFLO 
superconductor, the momentum of the Cooper pair in 
WSMs arises from the finite momentum of the WP and 
does not affect the energies of two paired electrons, b) 
The Weyl fermions are linearly dispersed and therefore 
their velocities are independent of their momenta. In or- 
der to guarantee the current conservation, only the con- 
tinuity of the wave function needs to be considered while 
the first derivative of the wave function does not help, 
so that momentum P^ carried by electrons and holes 
plays no role in the AR spectrum, c) The FFLO pairing 
states are formed within one Weyl node and thus the AR 
spectrum should also be independent of a, which is de- 
termined by a line connecting two WPs. Since the energy 
excitation is fully gapped and neither P± nor a makes 
contribution, the AR spectrum of an FFLO pairing WSM 
is isotropic and resembles that of a BCS pairing graphene 



The normalized AR spectrum for the BCS pairing 
WSMs is also plotted in Fig. [3] for comparison, which 
is a constant of 1.726 below the anisotropic effective gap 
| A sin a\ . It is clearly that the BCS and the FFLO pairing 
states are effectively distinguished by their AR spectra. 

It is helpful to discuss the experimental realization of 
the NS junction in the WSMs. The present effective 
model of the WSM is actually corresponding to the topo- 
logical insulator multilayers [I] and magnetically doped 
Bi2Se3 [5]. Several recent experimental progress on topo- 
logical insulators, such as the control of the mass term 
in BiTl(Si_5Se^) 2 [21) . the simultaneous magnetic and 
charge doping [22^, and the growth of the ultrathin high- 
quality films of Bi2Se3 [23] . paves a way to the realization 
of the WSM in topological insulator systems. The recent 
realization of superconductivity in the doped topologi- 
cal insulators [H] also indicates a hopeful future for the 
superconducting WSMs. 

In summary, the specular AR in the NS junction 
in inversion-symmetric WSMs is investigated. The 
AR spectrum for the BCS pairing states indicates an 
anisotropic effective gap |Asinct|, below which the nor- 
malized conductance is a constant of 1.726. The AR 
spectrum for the FFLO pairing states is independent of 
the momentum of the Cooper pair and is isotropic. The 
BCS and FFLO pairing states can be easily distinguished 
by their AR spectra. 
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